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We review quark model calculations of the transverse momentum dependent 
parton distributions (TMDs). For the T-even TMDs, we discuss the physi- 
cal origin of model relations which hold in a large class of quark models. For 
the T-odd TMDs we review results in a light-cone constituent quark model 
(CQM) with the final state interaction effects generated via single-gluon ex- 
change mechanism. As phcnomcnological application, we show the good agree- 
ment between results in the light-cone CQM and available experimental data 
for the Collins asymmetry. 
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1. Introduction 

Parton distributions entering many hard and exclusive processes play a 
key role to describe the nonperturbative structure of hadrons. The most 
prominent examples are the ordinary parton distributions which describe 
the probability to find in a fast-moving hadron a parton with a certain 
fraction of the longitudinal momentum of the parent hadron. However, in 
addition to the information on the longitudinal momentum distribution, a 
complete three-dimensional picture of the nucleon also requires knowledge 
of the transverse motion of partons. A full account of the orbital motion, 
which is also an important issue to understand the spin structure of the 
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nucleon, can be given in terms of transverse-momentum dependent parton 
distributions (TMDs). There are eight leading- twist TMDs. Two of them, 
the Boer-Mulders and Sivers functions are T-odd, i.e. they change sign un- 
der naive time reversal, which is defined as usual time reversal, but without 
interchange of initial and final states. The other six leading-twist TMDs 
are T-even. The latter will be reviewed in Sect. 2, discussing in particular 
the origin of model relations which were found in a large class of quark 
models 1-7 . Understanding the physical origin of these relations provides 
useful insights for building up realistic quark models. Indeed, from phc- 
nomenological applications to obscrvables one can verify to which extent 
the assumptions at the basis of the model relations are also satisfied in "na- 
ture". The T-odd TMDs will be the subject of Sect. 3, where we will show 
predictions in a light-cone constituent quark model (CQM) 8 which was 
successfully applied to describe many nonperturbative properties of the nu- 
cleon, like generalized parton distributions 9 , form factors 10 and distribution 
amplitudes 11 . The same light-cone CQM was also applied to describe single- 
spin asymmetries in semi-inclusive deep inelastic (SIDIS) which allow one 
to extract information on TMDs 12 . In Sect. 4 we will show as an example 
the results for the Collins asymmetry. 



2. T-even TMDs and quark model relations 

In the light-cone gauge A + = 0, the physical meaning of the correlations 
encoded in the T-even TMDs becomes especially transparent when using 
the canonical expansion in terms of light-cone Fock operators for the quark 
fields. Restricting ourselves to the quark contribution, one finds 

ft{x,k 2 x ) = (PA\V*(k)\PA), (1) 
Agl L (x,k 2 x ) = (PA\A*(k)\PA), (2) 

*±j£±gl T (x,k 2 ± ) = (PS ± \A«(k)\PS x ), (3) 

A^ h fl( x ,kl) = (PA\Tl(k)\PA), (4) 

S ± hZ T {xy ± ) + *±j^^ht2(x,ki) = {PS ± \T* x (k)\PS±). (5) 

where A (<Sj_) is the longitudinal (transverse) light-cone polarization of 
the nucleon, and P is the nucleon four-momentum. In Eqs. (l)-(5), the 
quark operators V q , A q , and can be written in terms of creation and 
annihilation operators of quarks with flavor q, light-cone helicity A and 
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light-cone momentum k = (xP + ,k±) as 



^sgn(A) 9 i(% A (fc), 



T?(k) 



(6) 



2ql(k)q_(k), 



where for a generic two-component vector a± we have defined clrx = a 1 ± 
ia 2 . The quark operators V and A q have a probabilistic interpretation since 
they are written just in terms of number operators N q = q\{k)q\(k). The 
operator T q ± has also a probabilistic interpretation but only when written 
in terms of transversely polarized operators 



9 =/2 s g ii ( s ±-HA k ) < is 



(k). 



(7) 



Choosing the frame system such that e z = P and e y = k±, the T-even 
TMDs can be classified in terms of the quark and nucleon light-cone polar- 
izations as summarized in Table 1. 



Table 1. T-even TMDs according to quark (columns) and nu- 
cleon (rows) light-cone polarizations. 
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Note: h^S = hi ± ijjh-h. 



In QCD all TMDs arc independent functions. However, in a large class 
of quark models 1-7 there appear relations among different TMDs. In par- 
ticular, for the T-even TMDs one finds i) three linear relations 



9\L 



2h\ 



0. 



91l 



9\t 



hi 



KZ = o, 



,-■-9 - 
1L 

2M 2 hlT 



with C q in Eq. (8) a real number, and ii) a quadratic relation 

2hif& = -(g q 1T ) 2 - 



(8) 
(9) 

(10) 



(11) 
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Although such model relations hold in quark models based on different 
descriptions of the quark dynamics, their origin can be traced back to a 
few common assumptions for treating the quark and nucleon spin degrees of 
freedom. Indeed, Eqs. (8)-(ll) connect the reciprocal polarizations of quark 
and nucleon, and we expect an underlying symmetry under rotation of the 
polarization degrees of freedom, common to all models. In order to describe 
the behavior of the polarizations under rotations, it is convenient to work in 
the basis of canonical (instant form) spin instead of light-cone helicity. This 
is because rotations in light-cone quantization depend non trivially on the 
interaction, while they become kinematical operators within instant-form 
quantization. A first common assumption to all the quark models satisfying 
the relations is that there is no mutual interaction among the quarks at 
the interaction vertex. In this case, light-cone helicity and canonical spin 
operators are just related by a Wigncr rotation in the polarization space 
with axis orthogonal to both e z and k± directions 

qx{k) = D% a > q s (k), 

where A and s refer to the light-cone helicity and canonical spin, respec- 
tively, while the explicit expression of the angle 9 in the different models 
can be found in Ref. [13]. Applying the rotation in Eq. (12) to the polariza- 
tion operators in Eq. (6) and the nucleon states, one finds the expressions 
of the T-even TMDs in the basis of canonical spin given in Table 2. If we 



Tabic 2. T-even TMDs according to quark (columns) and nucleon (rows) canonical po- 
larizations. 
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assume rotational symmetry around one of the coordinate axis, we find: 

• for the e x axis: the probability to find a quark with spin parallel 
to the longitudinal direction in a longitudinally polarized nucleon 
is equal to the probability to find both quark and nucleon trans- 
versely polarized in the e y direction. As a result, one can derive 
from Table 2 the two linear relations of Eqs. (9) and (10); 



n (i/ 2 )* _( cosf k L sin | \ 
Xa ~ UfcflSinf cosf I ' (iZ) 
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• for e z axis: the momentum distribution is axially symmetric when 
the canonical spins of the active quark and nucleon are aligned in 
the same transverse direction. Using the expressions in Table 2, one 
finds the quadratic relation of Eq. (11). 

The cylindrical symmetry around the e y leads to relations among TMDs 
which are not independent from the previous ones. The last linear relation 
in Eq. (8), relating unpolarized and polarized spin amplitudes, holds within 
more restrictive conditions. It requires spherical symmetry and 5*1/(6) sym- 
metry for the spin-isospin dependence. We note that all the models dis- 
cussed in the literature which satisfy the relations (9) and (10) among 
polarized TMDs and the quadratic relation (11) satisfy spherical symmetry 
in the instant-form, which is equivalent to assuming cylindrical symmetry 
around all the three spatial axis. Note that this condition is sufficient but 
not necessary for the validity of the relations. 

3. T-odd TMDs 

The gauge-link operator entering the quark correlation function which de- 
fines the TMDs is crucial to obtain non-zero T-odd quark distributions. In 
the light-cone CQM 8 , we consider the expansion of the gauge-link operator 
at leading order, taking into account the one-gluon exchange mechanism 
between the struck quark and the nucleon spectators described by (real) 
light-cone wave functions (LCWFs). This approach is complementary to a 
recent work 8 where the rescattering effects are incorporated in augmented 
LCWFs containing an imaginary phase. The approximation of truncating 
the expansion of the gauge-link operator at leading order is common to 
most of the quark-model calculations in literature 15,16 . However, studies 
beyond the one-gluon exchange approximation where one resums all order 
contributions have recently been presented in Ref. [17]. Following Rcf. [8], 
one can represent the T-odd TMDs in terms of overlap of light-cone am- 
plitudes corresponding to different orbital angular momentum. Both T-odd 
TMDs are obtained from the interference of wave components which differ 
by one unit of orbital angular momentum. In particular, we found that the 
Sivers function for both up and down quarks is dominated by the interfer- 
ence of S and P wave components, while the P — D wave interference terms 
contribute at most by 20%. On the other side, the relative weight of the 
P — D wave interference terms increases in the case of the Boer-Mulders 
function, in particular for the down-quark component. In Fig. 1 we show 
the results for the first transverse-momentum moments of the Sivers and 
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Fig. 1. Results for the first transverse-momentum moment of the Sivers and Boer Mul- 
ders functions for up and down quarks, as function of x. See text for the explanation of 
the different curves. 



Boer-Mulders functions. The dashed curves correspond to the results at 
the hadronic scale of the model. The solid curves are the results evolved to 
Q 2 = 2.5 GeV 2 . Since the exact evolution equations for the T-odd quark 
distributions arc still under study, we used those evolution equations which 
seem most promising to simulate the correct evolution. We evolved the first 
transverse-momentum moment of the Sivers function by means of the evo- 
lution pattern of the unpolarized parton distribution, while for the first 
transverse-momentum moment of the Boer- Mulders we used the evolution 
pattern of the transversity. After evolution, the model results are consistent 
with the available parametrizations. For the Sivers function, the lighter and 
darker shaded areas are the uncertainty bands due to the statistical error of 
the parametrizations of Rcf. [18] and Ref. [19] , respectively, which refer to an 
average scale of Q 2 = 2.5 GeV 2 . For the Boer-Mulders function, the dashed- 
dotted curves are the results of the phcnomcnological parametrization of 
Rcfs. [20,21] at the average scale of Q 2 = 2.4 GeV 2 , and the short-dashed 
curves show the results of Refs. [22,23] valid at Q 2 m 1 GeV 2 , with the 
shaded area describing the variation ranges allowed by positivity bounds. 

We also note that a non trivial constraint in model calculations of the 
Sivers function is given by the Burkardt sum rule 24 , which corresponds 
to require that the net (summed over all partons) transverse momentum 
due to final-state interaction is zero. Our model calculation of the Sivers 
function reproduces exactly this sum rule. 



4. Results for the Collins asymmetry 

In Ref. [12] the present results for the T-even TMDs were applied to esti- 
mate azimuthal asymmetries in SIDIS, discussing the range of applicability 
of the model, especially with regard to the scale dependence of the ob- 
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Fig. 2. The single-spin asymmetry A S j^^ h+ ^ s ^ = — A^j? c in DIS production of 
charged pions off proton and deuterium targets, as function of x. The solid curves are 
obtained using the light-cone CQM predictions for hi(x, Q 2 ) from Refs. [3,30]. The pro- 
ton target data are from HERMES 28 , the deuterium target data are from COMPASS 29 



scrvables and the transverse-momentum dependence of the distributions. 
Here we review the results for the Collins asymmetry A s ^^ + ^ s \ due to 
the Collins fragmentation function and to the chirally-odd TMDs hi. For 
the Collins function we use the results extracted in Rcf. [25]. In the de- 
nominator of the asymmetry we take fi from Rcf. [26] and the unpolarizcd 
fragmentation function from Rcf. [27], both at the scale Q 2 = 2.5 GeV 2 . 
In Fig. 2 the results for the Collins asymmetry in DIS production of charged 
pions off proton and deuterium targets are shown as function of x. The 
model results for hi evolved from the low hadronic scale of the model to 
Q 2 = 2.5 GeV 2 ideally describe the HERMES data 28 for a proton target 
(panels (a) and (b)). This is in line with the good agreement between our 
model predictions 30 and the phenomenological extraction of the transversity 
and the tensor charges. 31 Our results are also compatible with the COM- 
PASS data 29 for a deuterium target (panels (c) and (d)) which extend down 
to much lower values of x. 
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